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ABSTRACT 

We use a suite of hydrodynamical simulations of galaxy mergers to compare star 
formation rate (SFR) and black hole accretion rate (BEAR) for galaxies before the 
interaction (‘stochastic’ phase), during the ‘merger’ proper, lasting ^ 0.2 — 0.3 Gyr, 
and in the ‘remnant’ phase. We calculate the bi-variate distribution of SFR and BHAR 
and define the regions in the SFR-BHAR plane that the three phases occupy. No strong 
correlation between BHAR and galaxy-wide SFR is found. A possible exception are 
galaxies with the highest SFR and the highest BHAR. We also bin the data in the 
same way used in several observational studies, by either measuring the mean SFR 
for AGN in different luminosity bins, or the mean BHAR for galaxies in bins of SFR. 
We find that the apparent contradiction or SFR versus BHAR for observed samples 
of AGN and star forming galaxies is actually caused by binning effects. The two types 
of samples use different projections of the full bi-variate distribution, and the full 
information would lead to unambiguous interpretation. We also find that a galaxy 
can be classified as AGN-dominated up to 1.5 Gyr after the merger-driven starburst 
took place. Our study is consistent with the suggestion that most low-lnminosity AGN 
hosts do not show morphological disturbances. 

Key -words: 

galaxies: active - galaxies: interactions - galaxies: nuclei 


1 INTRODUCTION 


The correlation (or lack thereof) between the black hole ac¬ 
cretion rate (BHAR) and the star formation rate (SFR) 
of their host galaxies has been the subject of numer- 


Rosario et al. 2012 

-.g ., I- . 1 - , .. 

: Harrison et al.l 

2012- Rosario et al.l 

2013bl:lNetzer et al.l 

2 OI 4 I: iHickox et al. 

|2014|, and references 


therein). One of the drivers behind these studies is the em¬ 
pirical correlation be tween BH mass and bu l ge mass in the 
local universe fe.g.. [M arconi fc Hunt! l2003l : iHaring fc Rbd 
I 2 OO 4 I : iGiiltekin et al.l[2009ll which suggests that the BH and 
stellar bulge have assembled in tandem (co-evolution). In 
the strictest view of co-evolution, to obtain a BH mass pro¬ 
portional to the bulge mass, BHAR should be proportional 
to SFR (or, at least, the SFR that builds-up the bulge). In- 
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deed, the cosmic total SFR and BHAR (i.e., the rate per unit 
comoving volume) seem to track each other at least to z 3 


jHeAmanet_alJ^ 200j^ Merlom_et_^ 20041 : ISilverman et al.l 


I 2 OO 8 I. I 2 OO 9 I: iMadau fc DickinsoDll201.ib 


When SFR and BHAR are compared source by source, 
the connection appears to be weak, unless only the SFR 
in the central region is tak en into account (<1 kpc, 
iDiamond-Stanic fc Riek3 l2012l . and references therein). A 
comprehensive work of this type by iRosario et al.l ll2012ll . 
shows that at low AGN luminosities, BHAR and SFR are 
uncorrelated, while at high AGN luminosities a significant 
correlation emerges. On the other hand, Mullaney et al. 
(2012) and Chen et al. (2013) study a sample of mass- 
selected and star-forming galaxies respectively and suggest 
that the average BHAR correlates well with the average 
SFR, once the shorter variability time scales o f BHAR with 
respe ct to SFR are taken into account (see also lHickox et all 
l2014h . 
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2 Volonteri et al. 


In a companion paper (Volonteri et al. 2015) we have 
used new simulations of galaxy mergers and investigated the 
temporal correlation between SFR and BEAR, and their 
respective variability. We found that BEAR and nuclear 
(< 100 pc) SFR are well correlated and vary on similar 
timescales. Eowever, we found that gal axy-wide (< 5 kpc) 
SFR, which is u s ed in st atistical studies (iRosario et al.ll2012l: 
Mullanev et al.l I 2 OI 2 I : IChen et al.l I 2 OI 3 I : Delvecchio et al l 


2015h . and BEAR, are typically temporally uncorrelated, 


and have different variability timescales, except during the 
short-lived merger proper (~ 0.2 — 0.3 Gyr). We here com¬ 
pare SFR and BEAR by referring to the luminosity resulting 
from these processes (Lsfr and Lagn respectively), and aim 
at providing a framework for comparison with observations. 

In this Letter we show how very detailed numerical sim¬ 
ulation of galaxy mergers confirm the recent observational 
finding, provided the simulation results are treated in the 
same way as observations. In particular we show that ac¬ 
counting for the different durations of the various stages of 
the mergers, and binning the results in a way similar to 
the observational methods, can bring the results very close 
to the observed correlations and reconciles seemingly con¬ 
flicting results. The results present ed here are based on our 
previous detailed work published in ICapelo et al.l ([201^ and 
Volonteri et al. (2015). We refer the reader to those papers 
for a complete description of the simulations, of the physical 
implementation and of the analysis. 


2 COMPARISON WITH OBSERVATIONS 

The comparison of the numerical simulations with ob- 
served correlation s has some caveats. As discussed by 
iNeistein fc Netz^ (|2014| L measurements of SFRs are based 
on indicators that last >100 Myr, while the AGN luminosity, 
hence the BEAR, is measured instantaneously. Because of 
this, all our results concerning SFR represent averages over 
100 Myr. The SFR is converted to SF luminosity, Lsfr, by 
assuming that one solar mass per year corresponds to 10^° 
solar luminosities. BEAR is averaged instead over 1 Myr 
(see Volonteri et al. 2015 for a case where BEAR and SFR 
are measured over the same time interval). We distinguish 
between SF dominated galaxies {Lsf > Lagn) and AGN 
dominated galaxies {Lsf < Lagn)- 

Dense gas in the host galaxy may cause some obscu¬ 
ration, especially during the merger phase. We estimated 
the column density for a 1:2 merger, in cylinders of radius 
r = 100 pc centred on the BEs, and found that it varies 
between ~ 10^^ (face-on) and ~ 10^® cm“^ (edge-on), with 
an increase by a factor of 2 during the merger proper. Using 
the WebPIMMS interface (with a photon index F = 1.4), 
the average (maximal, for edge-on galaxies) dimming in the 
2-10 keV band varies between 20% (40%) in the stochastic 
and remnant phase, to 40% (55%) during the merger if the 
source were at z = 0, and it is almost negligible for 2 = 3. 
We discuss this further in section 2.2. 

An additional difference between theory and observa¬ 
tions is related to the fact that observed samples do not con¬ 
tain information about the fraction of sources in the stochas¬ 
tic, merger and remnant stages, and how many represent a 
merger with mass ratio of 1:2, 1:6 or 1:10. As shown in Fig.[T] 
the three stages cover differently the SFR-BEAR plane and 


a proper comparison with the observations requires an ad¬ 
justment of the various numbers in the three groups. 

The approach we adopt is to specify the boundaries of a 
region in the BEAR-SFR plane relevant to a population of 
relatively low-mass galaxies that includes quiescent galaxies 
and mergers with mass ratio from 1:1 to 1:10. Real cor¬ 
relations depend on the distribution of points within these 
boundaries, and will include galaxies of different masses and 
structural parameters, and this we cannot do quantitatively. 
One can, however, estimate the fraction of objects in each 
of the merger phases that is included in most observed sam¬ 
ples, e.g., the remnant phase in low-redshift samples, and 
remnants-l-mergers in higher-redshift samples. The duration 
of each phase must also be taken into account, since this de¬ 
termines the number of objects from this phase in a specific 
sample and hence the nature of the derived correlations that 
depend on such numbers. For example, the merger phase, 
being short-lived, ~ 0.2 — 0.3 Gyr, is sub-dominant in the 
general distribution, with respect to the stochastic and rem¬ 
nant phases, which last 0.8 — 1.3 Gyr each. 


2.1 Bi-variate distribution 


In Fig. [T]we show the bi-variate distribution of BEAR and 
SFR within 100 pc, and within 5 kpc, distinguishing galax¬ 
ies in the stochastic (red), merger (gold) and remnant (dark 
gray) phases. In the case of the central SFR, SFRioopc (left 
panel), the stochastic and remnant phases occupy similar re¬ 
gions, while the merger phase, beyond a primary peak coinci¬ 
dent with the other two phases, presents a secondary peak at 
higher BEAR and SFRioopc. The secondary peak is limited 
in its extent in SFR because averaging over 100 Myr removes 
the highest peaks of SFR, which are associated to the high¬ 
est peaks in BEA R (cf. Fig. 14 in Volon teri e t al. 2015). As 
noted previously (|Ej>pkins_J^Aluataeilj 2010|; ThnckCTet_^jJ 


2014 Diamond-Stanic fc Rieka 120121: iLaMassa et al.ll2oia 


Esauei et al.l 2014l ~) SFRioopc is a better tracer of BEAR 


than SFRskpc, i.e., they are be tter correlated. As discusse d 
by Capelo et al. 2015 (see also lEopkins fc Quataer3l201Cll L 
the main drivers of the BEAR are t he gas content (set¬ 
ting an overall ‘n ormalisation’, see also lRosario et aLll2013al : 
IVito et al]|2014h and local losses of angular momentum, and 
these became more clearly linked to SFR in the nuclear re¬ 
gion. We have compared the results of our standard run to 
a low-resolution run (see Appendix of Volonteri et al. 2015). 
The degree of correlation is degraded at lower resolution, 
showing that our simulations are accurately representing in¬ 
flows due to angular momentum loss. Additionally, because 
the nuclear region contains both cold inflowing gas and gas 
affected by thermal feedback, we calculate the accretion for 
each individual particle separately, rather than averaging 
over the gas broadly. The hot gas contribution is fairly neg¬ 
ligible, because of its lower density and higher temperature. 

The link between galaxy-wide SFR (SFRskpc, right) and 
BEAR seems much weaker. Galaxies in the stochastic phase 
occupy a region in the BEAR-SFRskpc plane (red contours 
in Fig. [TJ roughly tracing the region between the relations 
suggested by Mullaney et al. 2012 and Chen et al. 2013 to 
characterise star-forming galaxies (dotted blue lines). The 
merging systems (go ld contours) approac h the upper part of 
the curves shown in iRosario et al. H2 qT 3), who suggest that 
the highest luminosity AGN are merger-driven, and they 
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Figure 1. BHAR versus SFR within 100 pc (left), and within 5 kpc (right). Contours are based on ten equally-spaced logarithmic levels, 
dropping the six lowest levels for clarity. The SFR is an average over the 100 Myr before the time-step used for the BHAR calculation. 
We distinguish stochastic (red), merger (gold) and remnant (dark gray) phases. We show curves from Rosario et al. (2012, AGN, green, 
see section 3.2 for details), Mullaney et al. 2012 (mass-selected galaxies, upper blue dotted line) and Chen et al. (2013, star-forming 
galaxies, lower blue dotted line). We also mark the line separating AGN and SF dominated regions (light blue dashed line). 


exhibit a tighter correlation between BHAR and SFRskpc. 
In fact, in the merger phase, losses of angular momentum 
are driven by global dynamics and BHAR and SFR become 
somewhat better correlated even on large scales. In the rem¬ 
nant phase (gray contours) a wide range of BHAR can be 
associated to a given SFRskpc- 

The remnant phase is, perhaps, the most interesting 
when comparing to AGN observations, as 80 percent of 
the AGN do not show any hint of a companion (e.g., 

I Rosario et al ][2^ and are ordinary (massive) star forming 
galaxies. For most of the remnant phase the simulated galax¬ 
ies do not show strong morphological disturbances, however, 
the BHAR is sufficiently high at times that the galaxies en¬ 
ter the AGN dominated region of the BHAR-SFRskpc dia¬ 
gram. Observationally, this has been a matter of some de¬ 
bate. Rosario et al. (2013) find a strong connection between 
AGN activity and SFR, i.e. most AGN hosts are on the 
main-sequence of star forming galaxies llElbaz et al.ll200'7l : 
iNoeske et al.ll2f)(ra ). It has bee n suggested that th ere is an 
excess of AGN in post-merger (lEllison et al.ll2013l) or post- 
starburst galaxies dwild et al ]|2010l) . and that galaxies host¬ 
ing moderate-luminosity AGN are transitioning from SF to 
quiescence (Schawinski et al. 2009). Our broad interpreta¬ 
tion is that the AGN is more likely to be observed in the 
(much longer) remnant phase, where it enters the AGN- 
dominated region. In fact, in our calculations (Fig 13 in 
Volonteri et al 2015), the AGN zigzags continuously in and 
out of this region. 

We can now estimate the fraction of galaxies in different 
stages of the merger that are likely to be found in large ob¬ 
served samples. Galaxies in the merger phase would be con¬ 
sidered AGN-dominated (Lagn > Tsfr) for about 40 per 
cent of that phase, galaxies in the remnant phase are AGN- 


dominated for 25 per cent of their phase, and galaxies in the 
stochastic phase 7 per cent. However, the merger phase is 
much shorter. Therefore the overall probability, defined as 
the ratio between the time when Lagn > Tsfr and the total 
simulation time (stochastic -|- merger -|- remnant), of find¬ 
ing a galaxy in the remnant phase in the AGN-dominated 
region is almost twice as large than for a merging galaxy 
(13 per cent versus 7 per cent respectively), and six times 
higher than for a galaxy in the stochastic phase. 

2.2 BHAR versus SFR using binned data 

Most observational papers cannot consider the full bi-variate 
distribution of SFR and BHAR. Instead, they focus on the 
mean SFR for AGN in different luminosity bins (i.e. stack¬ 
ing in bins of BHAR; e.g., Rosario et al. 2012) or the mean 
BHAR in SFR bins (i.e. stacking in bins of SFR; e.g., Mul¬ 
laney et al. 2012, Chen et al. 2013 and Delvecchio et al. 
2015). To compare with those studies, we binned our simu¬ 
lations in the same way as used by them: i.e., we calculate 
the mean SFR in bins of BHAR and the mean BHAR in bins 
of SFR. Our points are not uniform on the SFR and BHAR 
axes across the different evolutionary stages, while obser¬ 
vational samples use relatively well distributed bins across 
their range. We therefore discard bins with few points (< 50 
for sub-samples, < 300 for the complete sample), and draw 
only 50 or 300 random points respectively for the remaining 
bins. 

The results are shown in Fig 2 where the first way of 
binning is shown in the top panel and the second in the 
bottom panel. Remarkably, in the former case we recover 
the trends of Rosario et al. (2012), while in the latter we 
recover the trend found by Mullaney et al. (2012), Chen 
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Figure 2. Mean SFR in bins of BHAR (top) and mean BHAR 
in bins of SFR (bottom); merger and remnant points have been 
slightly shifted to enable appreciating the size of the one-cr error 
bar. For the comparison, note that Rosario et al. (2012; green) 
measured the mean SFR for AGN in different luminosity bins, 
while Mullaney et al. (2012; blue) Chen et al. (2013; blue) mea¬ 
sured the mean BHAR, by stacking galaxies in bins of SFR. 

et al. (2013) and Delvecchio et al. (2015). We suggest that 
the different trends found for AGN and SF galaxies are in 
part caused by the different projections of the full bi-variate 
distribution, and that the intrinsic distribution of properties 
in those samples is similar to the one shown in Fig. [T] 

In light blue we show a bona fide global sample which 
takes into account the time spent in each stage, and where 
we have enforced high statistical significance in each bin. 
This is the sample to be compared to observations (SFR > 
0.5 Mq yr“^ seems also a plausible lower limit for the SFR). 
The process of averaging, and the statistical significance, 
have an obvious effect in removing the highest BHAR/SFR 
sources. For instance, when we require a large number of 
points per bin, the less populated high BHAR (SFR) bins 
disappear. Also, the mean in a BHAR (SFR) bin is driven 
by the more numerous sources at low SFR (BHAR). 

In the top panel, except for the two bins at the high¬ 
est BHAR, the SFR is consistent with being uncorrelated 
with BHAR. The wiggles show that the weight is dominated 
by the stochastic phase at low BHAR, then the remnant, 
and finally the merger phase take over. The bend at high 
BHAR is driven only by the last two points, dominated by 
the merger stage, as proposed by Rosario et al. 2012. How¬ 
ever, the increase in SFR is only a fact o r 2 to 3, similar to 
the enhancement seen bv IStanlev et al.l (l2015l ). An increase 
in the average masses of hosts of the most luminous AGN, 
or an enforced correlation between SFR and the long term 
BHAR may also be responsible for the trend. 

Note, again, that our simulations apply only to small- 
to-medium stellar mass systems. We can try to extrapolate 
the behaviour to more massive galaxies using the arguments 
outlined in Volonteri et al. 2015, that are based on the as¬ 


Figure 3. Trends of SFR and BHAR, for binned data, dividing 
the sample in AGN-dominated (darker points) and SF-dominated 
(lighter points). The BHAR at a given SFR for AGN-dominated 
galaxies is about 20 times higher than for SF-dominated ones. 


sumption that both BHAR and SFR would increase approxi¬ 
mately linearly with stellar mass. In this case, larger galaxies 
would move by about one order of magnitude of BHAR per 
one order of magnitude of SFR in the plane shown in Fig. [T] 
This leads us to suggest that if our stellar masses and BHs 
were two orders of magnitude larger (a few times ~ 10^^ 
and 10® solar masses respectively) the n they would occupy 
a region surrounding the 2 = 2 curve of iRosario et al.l (1201211 
when we bin in BHAR. In Fig. [2] we show the 2 = 0 curve 
shifted by 0.5 dex in both BHAR and SFR (the shapes of all 
these curves in Rosario et al. 2012 are basically identical). 

When binning in SFR (lower panel of Fig. (2)1, galax¬ 
ies in the merger and remnant phases are ~ 0.5 — 1 dex 
above the extrapolation of the relation suggested by Mul¬ 
laney et al. (2012) and Ghen et al (2013). The stochas¬ 
tic phase is closer to their fit in normalization, but has a 
shallower slope, caused by a drop in BHAR at the highest 
SFR, when, during a SF burst, SN feedback becomes strong 
enough to affect the gas properties. The full population has 
a best-fit log(BHAR) = (1.05±0.29) log(SFR)-2.27±0.12. 
If we included obscuration due to gas in the galaxy, the nor¬ 
malization decreases to —2.43. Furthermore, SF-dominated 
galaxies {Lsf > Lagn, light squares in Fig. [3]l have 
log(BHAR) = (0.93 ± 0.17) log(SFR) - 2.89 ± 0.07, slightly 
shallower and highly normalized than, but very close to, the 
curve proposed by Mullaney et al. (2012). Obscuration in 
the host would further decrease the normalization to —2.96. 

Beyond binni n g, sele ction also matters. As highlighted 
by Stanley et akl (l2015l ). their study and iRosario et al.l 
mm start with a population of AGN. Mass- or SF-selected 
samples start from a parent population of galaxies, which 
may or not host an AGN, and the average is likely dominated 
by sources which do not have AGN. Starting from an AGN 
population one studies the direct link between BHAR and 
SFR for galaxies in an active phase. A galaxy parent popu- 
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lation instead highlights the broad BHAR-SFR co-evolution 
over cosmic time and the whole population. In Fig. |3]we pro¬ 
pose a new way of looking at the SFR-BHAR connection. 
Observations often have only upper limits for either SFR or 
BHAR. In simulations, instead, we know when the BH is ac¬ 
creting or not, and what the SFR is, down to very low values 
(e.g., when SF is suppressed by SN feedback after a burst). 
We divide our sample in AGN- and SF-dominated galaxies 
{Lsf larger or smaller than Lagn respectively) and esti¬ 
mate the increase in BHAR versus SFR in galaxies during 
the phase where BHs are actively growing. AGN-dominated 
galaxies have a BHAR-SFR relation with a normalization 
about 20 times higher than SF-dominated ones, i.e. their 
BHs grow relatively more than the galaxy stellar mass by 
this factor. As the BHAR, i.e., the AGN luminosity limit, is 
lowered, more of the SF-dominated region is sampled, and 
at low luminosity AGN-selected samples cross from AGN- 
dominated to SF-dominated, thus lessening, or erasing, the 
significance of a correlation between SFR and BHAR (see 
also Hickox et al. 2014). 


3 CONCLUSIONS 

We analyse a suite of high-resolution galaxy merger sim¬ 
ulations to study BH and galaxy properties during various 
phases of the mergers. We calculate the BHAR and the SFR 
during the ‘stochastic’ phase (galaxies in isolation or in the 
early phases of an encounter), the ‘merger’ proper (when 
the merger dynamics dominates), and the ‘remnant’ phase 
(from the end of the merger to the return to quiescence). We 
find that in the remnant phase the BHAR can be sufficiently 
high at times to move the galaxy into the AGN dominated 
region. The probability of finding a galaxy in the remnant 
phase, even long after the starburst took place, in the AGN- 
dominated region is twice as large than for a merging galaxy. 
The main goal in this Letter is to compare the relationships 
between BHAR and SFR with observations. 

We find that different projections of the bi-variate dis¬ 
tribution recover different trends of the population. If the ob¬ 
servations can be extrapolated to the simulated mass range, 
we are able to reconcile seemingly contradictory observa¬ 
tional results when the mean SFR for AGN in different lu¬ 
minosity bins (e.g., Rosario et al. 2012) or the mean BHAR 
stacking galaxies in bins of SFR are measured (e.g., Mul- 
laney et al. 2012, Ghen et al. 2013, Delvecchio et al. 2015). 
Hickox et al. (2014), with an ansatz that over long timescales 
BHAR and SFR are perfectly correlated, reach similar con¬ 
clusions, also suggesting that future progress requires a di¬ 
rect measurement of the bivariate SFR/BHAR distribution. 

The bi-variate distribution derived from our simulations 
suggests that galaxies in the stochastic phase would not 
be considered AGN-dominated {Lsf 3> Lagn)- Galaxies 
would be AGN-dominated chiefly during the merger and 
remnant phases. The BHAR at a given SFR for AGN- 
dominated galaxies is about 20 times higher than for SF- 
dominated ones. During most of the merger proper and in 
the remnant phase a given SFR can be associated with a 
large range of BHAR. A possible exception is the groups 
of galaxies with the highest SFR and the highest BHAR, 
characterised by having Lsf — Lagn- 
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